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Abstract 



o 



Feldman and Karlin conjectured that the number of isolated fixed points for deterministic models of viability selection 
and recombination among n possible haplotypes has an upper bound of 2" — 1. Here a proof is provided. The upper 
bound of 3"~^ obtained by Lyubich et al. (2001) using Bezout's Theorem (1779) is reduced here to 2" through a change 
of representation that reduces the third-order polynomials to second order. A further reduction to 2" — 1 is obtained 
using the homogeneous representation of the system, which yields always one solution 'at infinity'. While the original 
conjecture was made for systems of viability selection and recombination, the results here generalize to viability selection 
with any arbitrary system of bi-parental transmission, which includes recombination and mutation as special cases. An 
example is constructed of a mutation-selection system that has 2" — 1 fixed points given any n, which shows that 2" — 1 
is the sharpest possible upper bound that can be found for the general space of selection and transmission coefficients. 

Keywords: Feldman Karlin conjecture; selection; recombination; transmission; fixed points; equilibria; Bezout's Theo- 
rem; homotopy method. 
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1. Introduction 



In a tribute issue to the late Sam Karlin, Feldman 



(20091 recounts their early collaborations (Feldman and 
KarUn[ |1968[ |KarUn and Feldman[ |1969[ |1970a|b[ ), and 



mentions a longstanding unsolved conjecture they proposed 
regarding the maximum number of isolated fixed points of 
the population genotype frequencies, under viability selec- 
tion and recombination: 

For the two-locus two-allele problem these con- 
siderations suggested a maximum of fifteen fixed 
points, and in our work with the symmetric 
viability model we demonstrated that fifteen 
was indeed realizable when recombination was 
present. Amazingly, to this day, our conjec- 
ture that the maximum number of equilibria 
in any rt-chroniosonie viability system and for 
any recombination arrangement is 2" — 1 has 
not been proven, although there are no coun- 
terexamples. 

Here I provide a proof, through a modification of the 
approach used by Lyubich (19921. The proof also gener- 



system, as described below. These results apply to sys- 
tems of autosomal loci with discrete, non-overlapping gen- 
erations, random mating, and constant viability selection. 

2. The Model 

The dynamical system considered here is represented 
by the recursion: 



x'r = 5j(x) := — 



1 



j,k=l 



TijkWjkXjXk (1) 



where 



iG{l,...,n} indexes the n possible gamete genotypes 
(i.e. haplotypes); 

xi, . . . , Xi, . . . Xn are the state variables, the frequencies 
of haplotypes in the population, so > and 

— Ij 



alizes the result to other genetic processes besides recom- 
bination, in fact, to any arbitrary biparental transmission 



\Xr, 



;(x) 



/.9i(x)^ 



\5n(x)y 
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x[ is the frequency of haplotype i in the next generation; 

Wjk = Wkj > is the fitness of the diploid genotype com- 
posed of haplotypes j and k; 

Tijk = Tikj > is the probabihty that diploid genotype jk 
produces gamete genotype i, so 

n 

Tijk = 1 for each j, k. 

i=l 

Perfect transmission is said to occur if gamete genotypes 
are identical to one or the other of the parental haplotypes 
in equal proportions: 



where Sij = 1 if i = j, dij = if i 7^ j. 
This recursion defines the map 

g : A""i ^ A"~i 

on the n — 1 dimensional simplex, 



A"-i = {x > 0: = 1} C M" 



T,, 



ijk 



of transmission 



The entries of the n by matrix 

probabilities are determined by the biological processes 
that occur during genetic transmission. Recombination, 
mutation, gene conversion, segregation distortion, inver- 
sions, and all of their combinations, are all representable 
by sets of Tijk- An absence of any transforming processes, 
and no segregation distortion, produces perfect transmis- 
sion. In this absence of transformation processes, a mul- 
tiple locus system is equivalent to a single locus system 
where each multilocus haplotype is formally a single al- 
lele. 

Processes that are not covered by this framework in- 
clude gene duplication and transposition, since the num- 
ber of potential haplotypes becomes infinite. Also, infinite 



alleles models (Ewens 2004, Sec. 3.6) are obviously not 



covered. Deletions, however, within a fixed set of genes, 
are representable by Tijk with finite n. 

The set of all possible transmission matrices clearly 
includes many that do not correspond to any known bio- 
logical processes, but this will be seen to be irrelevant, as 
the results apply to all transmission matrices. 

3. Prior Work 

Several previous studies prove results close to Feldman 
and Karlin's conjecture. Moran (1963), at the tail end of 



his Conclusions, provides early insight into the question of 
how many isolated fixed points are possible in a system 
of selection and recombination, when he invokes Bezout's 



The above discussion also raises the interest- 
ing theoretical question of how many station- 
ary points the adaptive topography for two un- 
linked loci can have. It is clear that in trivial 
cases such stationary points can fill up a lin- 
ear or a real continuum. This occurs when the 
Wij are independent of i, or j, or both. We 
may, however, ask how many isolated station- 
ary points are possible. The two equations typ- 
ified by (15) are cubics in P and p and hence, 
by Bezout's theorem have at most 9 distinct 
isolated solutions, real or complex. 



Tallis ( 1966 1 shows that with selection and perfect 
transmission, there are a maximum of 2" — 1 isolated fixed 
points of g, each fixed point being a population contain- 
ing haplotypes from one of the 2" — 1 possible subsets of 
n haplotypes, the empty set excluded. (Also see Lyubich 
Theorem 9.1.1.) 



1992 



Lyubich ( 1992 1 shows that with arbitrary transmission 
but with no selection, 2""^ is an upper bound on the num- 
ber of isolated fixed points of g. Lyubich et al. ( 2001 ) show 
that with arbitrary transmission and selection, 3"~^ is an 
upper bound on the number of isolated fixed points of g. 



Theorem (1779a I 



Lyubich 


( 


19921 and 


Lyubich 



et al. (2001 1 to obtain the upper bounds relies on Bezout's 



As origuially stated (translated from French): 



Theorem 1. \Bezout\ (ITT^ 47. p. 24) The degree of the 
final equation resulting from an arbitrary number of com- 
plete equations containing the same number of unknowns 
and with arbitrary degrees is equal to the product of the 
exponents of the degrees of these equations. 

Here, the 'final equation' is a reference to the elim- 
ination method for solving systems of polynomials. An 
immediate application is: 



Corollary 1. Bezout (1779b. 48.3. p. 24) The number 



of intersection points of three surfaces expressed by alge- 
braic equations is not greater than the product of the three 
exponents of the degrees of these equations. 

The corollary is readily generalized to the intersection 
of arbitrary numbers of algebraic curves, to yield the ver- 
sion of Bezout's theorem utilized here, restated fromjKoMr 
(2008| p. 365-366): 



Theorem 2 (Bezout's Theorem). 

Let /i(x), . . . , fn{^) be n polynomials in n variables, and 
for each i let di be the degree of fi. Then either 

1. the equation(s) /i(x) = ••. = fn{x) = have at 
most did2 ■ ■ • dn solutions; or 

2. the fi vanish identically on an algebraic curve C , and 
so there is a continuous family of solutions. 



Theorem (1779a I: 
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3.1. In the Absence of Selection 

Bezout's Theorem is applied as follows by |Lyubich| 



(1992 pp. 294-296). In the absence of selection, the sys- 
tem of fixed points of ([ij can be written as the zeros of a 
system of polynomials: 



= /i(x) := ^ TijkXjXk ~Xi {i^ 1, 

3.k=l 



and 



= /„(x) :=^x, -1. 



1) (2) 



(3) 



i=l 



This gives n~l equations of degree 2, and one equation of 
degree 1 ( ([2]) and (|3|, respectively). Therefore, f(x) = 
can have at most 2" * 1 = 2"~^ isolated solutions, which 
are the fixed points of g. 



Lyubich (1992 Theorem 8.1.4 and Corollary 8.1.7, pp. 



295-296) invokes the Poincare-Hopf index theorem to re- 
duce the upper bound to 2"^^ — 1 when g maps all points 
on the boundary in the direction of the interior of the sim- 
plex: i.e., 

for x e aA"-! : 

(1 - e) X + e g(x) e A"-i - dA''-^ for small e > 0, 

or, more simply, 

if = 0, then ^^(x) > 0. (4) 



Lyubich ( 1992 , p. 295) provides inferences between several 
related conditions. Let ^-"(g) = {x: g(x) = x, x e A"^^} 
represent the set of fixed points of g on A"^^. The condi- 



tions are: 



^(g)n9A"-i =0; 

1 = Poincare-Hopf index of g on 9A"~^; 



T 



is an irreducible matrix; 



(5) 
(6) 

(7) 

(8) 
(9) 



Tijj > V i, j e {1, . . . 

g(x) e A"-i - aA"-\ Vx e A"-^. 

Lyubich proposes the following chain of implications: 

(|9|^(|8|=^(|7])=^([5])=>(|4])=4>([6]| 

Additional consideration, however, shows the correct 
chain of implications to be: 



(© ^ (§ 







In particular, ^ ([4|. Lyubich writes, if "the bound- 

ary of the simplex contains no fixed points [([5])] then the 
vector field Vx — x cm the boundary is directed inside the 



simplex [(|4])]." Following is a (nonbiological) counterex- 
ample where ([5| and (|7| hold, but ([4| does not (without 
claiming the theorems themselves to be invalid): it is pos- 
sible for the boundary to be without fixed points, and for 
the matrix \Tijj\ to be irreducible (Lyubich uses the alter- 
nate term 'indecomposable'), yet for g to map boundary 
points toward the boundary and not the interior. 

Let T simply rotate the index of each haplotype by -1: 



T,, 



ijk 



1 for i (mod n) 

1/2 for i (mod n) 

1/2 for i (mod n) 

otherwise. 



An initial point x = ei := (1, 0, • • • , 0)^ (t, the transpose) 
is taken by g through a cycle of period n through the 
vertices of the simplex, ei, e„, e„_i, . . ., 62, ei, e„, . . .. 

Every point on the boundary maps to a different point 
on the boundary. This can be seen because every boundary 
point must have some indices + l (modulo n) such that 
Xi > while Xi+i = 0, but T rotates the indices so that 
gi(x) = in the next generation, so no boundary point 
is fixed. Moreover, for any point on the boundary with 
adjacent zeroes, i.e. Xi = Xi^i = 0, then (/^(x) = 0, con- 
trary to boundary condition Here, boundary points of 
a sub-simplex map to other boundary points of that sub- 
simplex, so (1 — e)x -t- eg(x) remains on the boundary for 
all < e < 1. Any fixed points must therefore be in the 
interior of the simplex, and by symmetry, this can only 
be the center, x = (1/n, . . . , 1/n)^. Thus, the exclusion 
of fixed points from the boundary does not imply that g 
maps the boundary in the direction of the interior. 

What Lyubich is really after, however, is condition 
So if one starts by assuming ^ , then ^ follows and the 
rest of Lyubich's proof goes through: 



1. 



g has a Poincare-Hopf index of 1 on the boundary of 
the simplex ([6|; 

The index must equal the sum of the indices of the 



fixed points in the interior of the simplex (see Glass 
|1975b| and |Glass]|1975a| for an accessible explication 
of this approach); 

3. The index of a non-degenerate fixed point must be 
+1 or -1; 

4. Supposing that the maximum of 2"^^ fixed points is 
attained, then each must have multiplicity of 1, and 
is thus non-degenerate; 

5. Having 2"~^ non-degenerate fixed points in the in- 
terior of the simplex, however, would produce a sum 
for their indices that is even, contrary to the index 
of g under condition ^ ; 

6. So there can be no more than 2"^^ — 1 isolated fixed 
points of g in the interior of A"^^ given ([4|. 

3.2. In the Presence of Selection 

The above result is derived when selection is absent, 
and the only force acting is transmission of some arbitrary 
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form. When selection is included along with transmission, 
the system of fixed points becomes: 



= /.(x) 



TijkWjkXjXk - Xi WjkXjXk (10) 

j,k=l j,k=l 
(i = l,...,7i- 1), 



and 



= /„(x) 



The terms Xi 5^2fc=L '^jkXjXk in ( 10 ) are all of degree 3, so 



Lyubich et al. (2001 eq. (33)) apply Bezout's Theorem to 
obtain an upper bound of 3"~^ on the number of isolated 
fixed points. 

The vastly larger value for the upper bound when selec- 
tion is present, 3"^^, versus 2"^^ when selection is absent, 
seems counterintuitive, and one suspects that 3"~^ can be 
sharpened. 

4. Results 



Closer examination of ( 10 ) finds that the equations are 



third order only due to the presence of a single common 
factor, the mean fitness w :— Y^" k=i^jkXjXk- Such a 
structure suggests potentials from a change in the repre- 
sentation. This is indeed the case, and the system can be 
made second order by introducing an additional variable 
y, and an additional equation that constrains y to equal 
the mean fitness at equilibrium. The result is as follows: 

Theorem 3 (Generalized Feldman-Karlin Conjecture). 

Consider the evolutionary system with viability selection, 
random mating, and arbitrary transmission of n possible 
haplotypes, represented by the map g; 

1 " 

wjkXjXk ^•''^^ 

where i = 1,. . . ,n, Tijk > 0, X^Li ^ufe = ^' ^ij - ^' "'^'^ 



> 0, = 1. 



(11) 



The number of isolated fixed points of g is never greater 
than 2" - 1. 



equations in n + 1 variables: 

n 

= /i(x, y) := ^ TijkWjkXjXk - x^y (12) 
(i = 1, . . . ,n- 1), 

n 

= /„(x,2;):-^.T,-l, (13) 

n 

= fn+i{^,y) ■■= ^ WijXiXj - y. (14) 
Since all fixed points are isolated by hypothesis, applica- 



tion of Bezout's Theorem to (12 1, (13 1, and (14 1 gives an 
upper bound of 2"^^ * 1 * 2 ~ 2" on the number of fixed 
points. 

This comes within 1 of the value conjectured by Feld- 
man and Karlin for the upper bound on the number of 
isolated fixed points of evolutionary systems ([T|. Their 
upper bound could be demonstrated if one could show that 



a solution satisfying ( 12 1 , ( 13 ) , and ( 14 1 always exists out- 



side the simplex. Complex-valued solutions, and solutions 
'at infinity', qualify and must be included in the count of 
solutions in Bezout's Theorem. 

A solution 'at infinity' is formally accounted for by 
mapping the system to projective space; this is accom- 
plished by a homogeneous form for the intersection equa- 



tions (Shafarevich 1994 pp. 16-21). A new set of variables 
is defined: 



z, 

Xi 

y ■ 



- Xi z, and 
yz. 



Substituting in (12 1, (13), and (14), one obtains n -|- 1 



equations in n + 2 unknowns: 



1 " 

0=,A(x,y) E 



1 



(i = 1, . . . ,n- 1), 
1 " 

= /n(x,y) := - ^Xj - 1, 



i=i 



1 " 

= /„+i(x,?/) := ^ X! ^ii'^i^o -y/z- 



(15) 

(16) 
(17) 



Multiplying both sides of (15) and (17) z , and (16) by z, 



Proof. With the introduction of an additional variable y, 
the system of fixed points can be represented as n -I- 1 
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one obtains the homogeneous form: 

n 

= fi{x.,y,z) := ^ TijkWjkXjXk -xty (18) 
(i = l,...,n- 1), 



= /„(x,y,z) := 



(19) 



= fn+i{^^y,z) ■■= ^ w,jx,xj -y z. (20) 

In the homogeneous representation, any non-trivial solu- 
tion, (x, y, z) 7^ (0,0,0), gives as solutions all its scalar 
multiples c(x, y, z), c G C. Hence all scalar multiples of 
a solution count as a single point (a point in projective 
space) when counting solutions. 



The original system (12 1 is obtained by setting z = 1. 
By setting z = 0, any nontrivial solution (x 7^ or y 7^ 0) 
corresponds to a solution 'at infinity' for (12), (13 1, and 



(14 1, since it would give x = x/0 and y = y/0. 

Letting z = and x = 0, we see that this is a solution: 

n 

7,(0, 0) - ^ T,,kW,kO * - Oy - 
(i = 1,. 1), 

n 

7„(0,y,0) =^0-0-0, 

n 

7„+i(0, y, 0) = ^ w;^ * * - y = 0. 

ij = l 

The variable y is clearly unconstrained here, hence, the 
non-trivial solution is: 

(x,y) = c(0,l), c^O. 

This solution 'at infinity' reduces by 1 the number of pos- 
sible finite-valued fixed points, leaving an upper bound of 
2" — 1 on the number of isolated fixed points in the simplex 
A"-i. □ 

An illustration of how the solution 'at infinity' arises is 
given in section [5. 1[ 



5. A Mutation-Selection System Bearing 2" — 1 
Isolated Fixed Points 

Here, I show that 2" — 1 is the smallest possible up- 
per bound on the number of isolated fixed points for the 
general space of selection-transmission systems ([T|) , by con- 
structing an example that attain this bound. It is already 
known that the bound is attained by examples of systems 
with selection and perfect transmission, where one fixed 
point is located in the interior of the simplex and in the 
interior of each sub-simplex on the boundary, including 



the vertices (Tallis 19661. So what remains to be deter- 
mined is whether this bound can also be attained under 
some form of imperfect transmission. 

Such systems can be produced through the homotopy 



continuation method ( Kotsffeas 2001 Li 2003 1, where one 



creates a continuous family of systems between a known 
system, and an unknown system with desired properties. 
In this case, the homotopy will go from a known system 
of selection and perfect transmission that has 2" — 1 fixed 
points, to systems with imperfect transmission, by per- 
turbing the transmission probabilities. Under proper con- 
ditions, the homotopy will produce paths from the fixed 
points of the known system to the fixed points of unknown 
systems. 

The homotopy continuation method per se originated 



Drexler 


(1977 


), and 


1 — 

Chow et al. 


(1978 


1 ( 


Li 


1997 



An es- 
sential part of this method can be found earlier in the 
'method of small parameters' of Karlin and McGregor 



(1972b 1972c 1972a I: 



Principle I [if a system of transformations act- 
ing on a certain set (in finite dimensional 
space) has a "stable" fixed point, then a slight 
perturbation of the system maintains a stable 
fixed point nearby] can be interpreted as a per- 
turbation or continuity theorem. Starting with 
a given genetic system for which the nature of 
the equilibria can be fully delineated (for exam- 
ple, the classical multi-allelic viability model), 
it is desired to investigate a perturbed version 
of the model. (|1972b[ p. 86) 



Karlin and McGregor (1972c Theorem 4.4, p. 231) use 



the implicit function theorem to establish the existence, 
uniqueness, and nearness of fixed points under perturbed 
models. The stability properties of the fixed points also 
are preserved by the additional assumption that all fixed 
points are hyperbolic. Hyperbolicity is not needed, how- 
ever — only non-degeneracy — to establish the existence, 
uniqueness, and nearness of isolated fixed points ( |Akin[ 



1983 



p. 24). 

For the choice of the system to perturb, I use ([ij 

with perfect transmission, which is formally a one-locus, 
multiple-allele system: 



x'^ = ft(x) 



w^jXj (21) 



1 fixed 



Tallis (19661 gives two examples of (211 with 2 
points, a fully overdominant and a fully underdominant 
system. The fully overdominant system is due to [Wright | 
( 19491 (discussed in [Li[ [1955| 



P- 



Si, Wij = 1 for all i ^ J 
underdominant system has w. 
i ^ j, where < < 1. 



260) 
where < 
1 



and has wa = 1 — 
s, < 1. The fully 

Wij = 1 — s for all 
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The homotopy is created by perturbing perfect trans- 
mission using uniform mutation: 

T^Jk - (1 - M) ^ +S^k)+ M/n, (22) 

where fi is the mutation rate. With /i > 0, (22l satisfies 
conditions Q and yielding conditions (4f^ (|5| and 
([6|. Thus none of the boundary fixed points under ( [2l] ) 
can remain on the boundary for /i > 0. 

To keep all 2" — 1 fixed points in the simplex when 



H > 0, the fixed points of (21) on the boundary need to 
move inside the simplex. Since g points into the simplex 
for /i > 0; toward these fixed points, they need to have 



stable manifolds ( jWiggins 19901 pp. 193-239) that enter 
the simplex. A situation that produces this is where all 
n corners of A"~^ are stable sink nodes, and all fc-allele 
(fc G {2, . . . , 71 — 1}) polymorphic boundary fixed points 
are unstable source nodes with respect to the sub-simplex 
A'^"^ for which they are interior points. Kingman (1961 
p. 578) showed that under (21 1, the vertex equilibria are 



stable and all polymorphic equilibria unstable if W is posi- 
tive definite. This is the case with the fully underdominant 



system (Also see Christiansen||1990 ) 



Let the fitnesses be Wi 



1 



s Sij, where s > 0. 



In 




(21 1, the fixed points comprise: the vertices of the simplex; 
the centers of each h — 1 dimensional sub-simplex on the 
boundary, and the center of the simplex. Fixed points will 
be of the form 



(23) 



in any of (^) distinct permutations of the order, where h 
varies from 1 to n. Hence the number of fixed points is 

ELi(p = 2"-i- . 

Letting the mutation rate ^ become positive, this sys- 
tem provides a constructive proof for the following: 

Theorem 4. The upper bound of 2" — I on the num- 
ber of isolated fixed points for evolutionary systems ([ij 
is sharp over the general space of selection and transmis- 
sion coefficients. In particular, for any n, a system with 
imperfect transmission can be constructed that attains the 
upper bound o/ 2" — 1 fixed points in the simplex, using 
Wij = 1 -\- s Sij , s > 0, and 



Tijk = (1 - + (5jfe) + /^/n, 



for the ranges < /i < 1/2 and s > 



I/m-2 



, yielding: 



Proof. For small /i, (24) is a perturbation of the system 



with perfect transmission, so its fixed points will be close to 



the isolated and non-degenerate points ( 23 1 . By symmetry 
we can expect them to be of the form: 



i = P( 



1 - a 



1 - a 



(25) 



for h = l,...,n— 1, where P is any permutation of the 
order of the entries. For h — n or h — 0, Xi — 1/n ioi all 
i. 

To verify this form, and solve for a, we substitute into 



(24) 



1-a 



(1 + sV) 



,1 — a n , ,/l — a 

(1-m)^— + - + 

h n \ h 



V 



where 



a;? = 



(l-«)^ 
{l-af 



+ {n-h) 



(n-h) 



(n - hy 



which yields three roots: 

a = 1 — h/n, 

and 



1 /I 

"^2+^ 2 



± 



1 2 



fj,{h + s){n — h) 



The first root yields the central equilibrium Xi = 1/n. Sub- 
stitution shows that the other two roots yield the two equi- 



librium forms in (25 1, where h and n — h are interchanged. 
The term inside the radical must be non-negative for real 
solutions, and so imposes constraints on and s: 



l{h) := 



^{h + s){n — h) 



ns 

fin fi 
4^ ~ 2 



(26) 



> 0. 



To obtain the conditions for jih) > 0, first the value 
of h is found that minimizes "f{h): 



5«(x) 



En 
.7 = 1 • 



(24) 



(1 - [^i)xi ^ h s(i — [i)Xi ^ — ^ 

n n 



d 

dh 







h 


n 




n 



H+ -] fi 



So jtj{h) = at h ~ n/2. The second derivative, 
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Figure 1: A plot of the polynomials l |12[ l for the example |24j with 
n = 2, s = 1, wii = W22 = 2, u>i2 = W21 = 1, and fi = 0.05. The 
range is truncated at to show the curves of solutions to fi{xi,y) = 
0. Shown are the three points of intersection, where J/'*^^) = 



f2{i:\''\y^''^) = 0, for k = 1,2,3. At those three points, w 



Counting as the fourth intersection 'at infinity', 'Xj*'' indicates where 
the two curves approach parallel lines as j/ — > — oo. 

is positive, hence is a minimum: 

2 



7(n/2) = 



1 1 

1 /in /t 

4 ~ 47 ^ 2 



1 fin fi 

4 ^ 4^ ~ 2 



Therefore, 



7(n/2) > 



1 / n 
4 - ' V 4s 



/i < 1/2 and s > 



fj,n 



If n is odd, then the requirement is 7(^^^) = 7(^^^) > 0, 
but since j(n/2) is the minimum, the above constraints on 
/i and s are sufficient to keep 7(^^) non-negative. In ( pG] ), 

since fj.{h+s){n-h)/{ns) > 0, then y/^{h) < - 
so a > 0, as required. 

The number of fixed points is the number of permu- 



tations of (25 1, 2" — 1, the same as for (23) with perfect 



transmission. 



Therefore, g in (24) is an example of a system with 



imperfect transmission that attains the upper bound of 
2" — 1 isolated fixed points, all in the interior of the sim- 
plex. Thus, for the general space of fitness and transmis- 
sion coefficients, the upper bound in Theorem [3] is sharp 
and cannot be lowered. □ 



5.1. The Solution 'At Infinity' 

This example system is also useful to understand the 
nature of the 2"-th solution 'at infinity' in Bezout's Theo- 
rem. Fig. [l] shows the polynomials (12) for the two-allele 
case. The two algebraic surfaces intersect each other and 
the fi ~ plane at 2^ — 1 = 3 points. The absence of 
a fourth point is because the two algebraic surfaces ap- 
proach parallel lines along the — plane as y goes to 
— oo. These parallel lines are said to 'meet at infinity' 
when the system is represented in projective space using 



the homogeneous form ( 20 1 



6. Discussion 

Theorem [3] proves the conjecture of Feldman and Kar- 
lin on the maximum number of isolated fixed points in a 
system of selection and recombination, and extends it to 
arbitrary transmission processes, of which recombination 
and mutation represent special cases. Thi s substantially 
sharpens the previous upper bound of 3"~^ (Lyubich et al. 



2001 1 on the number of isolated fixed points of an evolu- 



tionary system with selection and arbitrary transmission. 

No attempt has been made to characterize the condi- 
tions on Tijk and Wjk that produce only isolated and non- 
degenerate fixed points. More on this issue can be found 
in Lyubich et al. (2001). One may mention, however, that 
such conditions are generic, in that for 'almost all' sets 
of n algebraic hypersurfaces of degree n, the intersection 
consists of isolated non-degenerate fixed points ( ,Shafare^ 



vich 



1994 



P- 



223, 



Garcia and Li 1980); and sets of Ti 



^ ijk 



values that produce degenerate fixed points are nowhere 
Theorem 



dense (Lyubich 1992 Theorem 8.1.3) in the space of Tijk 
values. Recombination-only systems without selection ex- 
hibit only continua of fixed points with no isolated fixed 
points, because arbitrary allele frequencies are all invari- 
ant, and only linkage disequilibria change in time. Certain 
non-generic selection regimes also yield continua of fixed 
points. 

The present paper does not touch at all upon the ques- 
tion of the stability of the equilibria. Studies on the 
stability of equilibria include the papers cited previously, 
and, as a selected additional list, Feldman et al. (1974 1, 
Karlin (1975), Karlin and Liberman ( 1976 1 , jFeldman and 



(19811, Hastings 



(19851, 



Puniyani and Feldman] ( 2006 ) 



Liberman| (|1979 l |Karliii|ll979J, Karlin] (|1980j, Hastings 



Franklin and Feldman ( 2000 1 , and 



It is notable that when selection is removed from the 
evolutionary system ([ij, the upper bound on the number 
of isolated fixed points decreases by half, from 2" — 1 to 
2^-1 (2"^i - 1 when condition (|4| holds). In contrast, 
removing arbitrarily complicated imperfect transmission 
from (jlj to constrain the systems to form (21 1, where only 
selection acts, does not reduce the potential number of 
isolated fixed points at all. 

Removing selection means setting the n (n + l)/2 in- 
dependent values of Wij to 1, whereas removing imperfect 
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transmission means setting the — 1) n (n + l)/2 indepen- 
dent values of T^-fc to {Sij + Sik)/2. Hence, paradoxically, 
when the 0{n^) fitness coefficients are allowed to vary, 
it doubles the potential number of isolated fixed points, 
whereas when the 0{n^) transmission probabilities are al- 
lowed to vary, the potential number of fixed points remains 
unchanged. This reveals a structural difference between 
the roles of selection and transmission in the dynamics. 

A homotopy of a continuous family of evolutionary sys- 
tems can be defined between any two evolutionary systems 
([ij, in particular between ones with and without selection 
through the parameterizations: 

Wij{a) = [1 — a) + a Wij, or Wij(a) — w°j. 

So any system 

1 



^ Wjk{a) Xj Xk ^■''-'^ 



Xj Xk 



with greater than 2"^^ fixed points will exhibit bifurca- 
tions of fixed points for some value(s) of a € [0, 1]. 

One can contemplate many other possible applications 
of the homotopy method to investigate the equilibrium 
behavior of these evolutionary systems. 
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